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Abstract 

In planar QCD, in two space time dimensions, the meson eigenvalue equation has a nonlocal 
structure interpretable as resulting from hidden degrees of freedom. The nonlocality can be recon¬ 
structed from the functional form of the pion mass dependence on quark mass within an expansion 
starting from a special one dimensional Schrodinger problem. The one dimensional problem makes 
the pion mass depend on the quark mass through a simple quadratic relation which is shown to be 
compatible also with numerical data obtained in four dimensions. 
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INTRODUCTION 


Intuitively, if strings indeed describe QCD, the best place to make this precise is in the 
meson sector, at inhnite number of colors, N. The mesons do not influence the gauge held 
vacuum and are open string probes of the closed string background describing it. The mesons 
are free then, and, for hxed quark mass (assume all quarks have equal mass for simplicity) 
their masses squared should fall into some regular patterns, extending to inhnite mass. An 
example is provided by SU{N) gauge theory at inhnite N in two space-time dimensions. 
Although in two dimensions the gauge degrees of freedom are unlikely to provide a rich 
enough structure for a fully featured closed string background, a degenerate form, closer to 


topological string theory, might exist 


Q. 


A speculation for four dimensions condenses all the unknown structure into one real 
function of one real variable, playing the role of a warp factor in a hve dimensional metric. 
This function determines the closed string background and the meson mass dependence on 


quark mass. In simple examples the warp factor directly determines the masses of some 
particles by a local partial diherential equation, but in the case of QCD it is unlikely that 
the equation giving the meson masses indeed is local in any set of variables. However, if 
one expands in energy, relative to a scale set by the string tension, the mass of the lightest 
mesons might be well described by a leading approximation, which does consist of a simple 
second order partial diherential equation. Higher orders in the two dimensional held theory 
describing the four dimensional QCD string will induce further corrections in the meson- 
quark mass dependence. As we shall show below, in two dimensions, the function giving 
the mass of the pion, m^, as a function of the mass of the quark, m^, starts from a simple 
quadratic formula, structurally rooted in an ordinary second order eigenvalue problem. The 
full, exact dependence then determines all the higher order corrections in a low energy 
expansion which reproduces the full non-local meson equation. This leads us to a numerical 
test of a similar quadratic approximation in four dimensions. 

The basic philosophy we adopt is to assume that mesonic spectral data could be used to 
reconstruct the unknown free string theory purportedly describing planar QCD, somewhat 
akin to an inverse scatteringapproach using spectral and scattering data to get the potential 


B. 


in a Schrodinger problem 

AdS/CFT motivated modeling seems to be relatively successful 


Q 


phenomenologically. 
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although for four dimensional planar QCD we only have ad-hoc motivated equations, and we 
even don’t yet have the numerical values of the masses with their dependence on quark mass 
we should compare these to, since the real world has only 3 colors. As already mentioned, 
it is unlikely that such equations, which work exactly in some AdS/CFT cases, at all exist 
as exact representations at infinite N, and it is much less likely that any such equations are 
exact for N = 3. 

In this letter we take some steps to improve our understanding of mesonic planar QCD 
in two and four dimensions by focusing on the dependence of the lightest pseudo-Goldstone 
mesons (pions) on the quark mass ai N = oo. 

Our work also provides a test of the numerical methodology we have been developing for 
dealing with planar QCD in the meson sector. Here, we shall skip most technical details, 
and concentrate on presenting the results, speculating on their possible meaning. 


TWO DIMENSIONS 


We start in two dimensions, from ’t Hooft’s exact solution. The chiral condensate (V’V’) 
and the meson masses are exactly known as functions of m^, the quark mass. The ’t 
Hooft coupling A = is used to set the scale. The physical dimensionless quantities are 
7 = A-i Trig and jj? = \ Spontaneous symmetry breaking occurs when the 

order of Urrrite is: leavirrg a .ro.r-.ero condeasafe ia the .ero quark mass 

limit p. 

The meson spectrum and its dependence on quark mass (here the quarks have been taken 


of equal mass for simplicity) is governed by’t Hooft’s equation j?]: 

fi (p{y) - 


7 
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dy = y (j){x) 


( 1 ) 


X is the fraction of meson light-cone momentum carried by one of the quarks and varies 
between 0 and 1. 

It is natural to ask whether ’t Hooft’s equation contains any structural hints that it might 
admit a geometrical interpretation tied to self-consisted string propagation, which also could 
employ extra dimensions for the string to propagate in. More precisely, as a first step, we 
wish to see whether the addition of some auxiliary continuous arguments can turn’t Hooft’s 
equation into a local differential equation of second order. The answer is positive. 
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While the hrst term in m is local in x, the second is highly non-local. One can easily 
expand it in derivatives of cp. To deal with the combinatorics of the coefficients in a more 
efficient way, we seek a set of functions that diagonalize this term. This set is easily found, 
using the integral: 


P dx 


{x — ay ^{b — xY 


7r(c — ay ^ 


■ cot(i/7r). 


( 2 ) 


Ja X — C {b — cY 

where a < c < b and 0 < < 1; we shall use the formula also at = 0,1. It is now 

convenient to introduce s = logj^, measuring the rapidity difference between the quark 
and anti-quark, and the equation gets recast into the following form: 

,2 
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( d\ 

1 — TT— cot 


ds 

\ dsj 


T(s) + 




4 cosh"^ 


-\k(s) = 7 \k(s) 


( 3 ) 


In this form of the equation, /i^ is more naturally viewed as part of the potential term, 
while the eigenvalue role is more naturally taken up by 7 , measuring the quark mass. Using 

^2 


Tip coth 7rp — 1 = 


p- 


^0 P + ^ 


( 4 ) 


it becomes obvious that we could localize the equation by adding a new dimension corre¬ 
sponding to the internal coordinate along and open string, where the modes are labeled by 
nonzero integers n. 

Explicitly, we introduce a real held x(s, a), where a G [0, tt], in addition to the held T(s), 
which plays the role of the meson wave function, but is taken as real. The held y obeys 
Neumann boundary conditions in a and is further constrained to contain no zero mode. 


%(5,0) ^ dx{s,7i) 
da da 


j X(s,a) = 0 


( 5 ) 


We now introduce a local action functional S'[T, y] whose extremization subject to the above 
boundary condition produces the meson wave equation. 


S['p, \] = la a) 


9^ I 9^ 
ds^ ' 




4 cosh^ (I) 


vl/(s) 


X(s, cr) = Zl cos(ncr)x„(s) Xn = X-r 

n^O 


( 6 ) 


Eliminating the variables Xn{s) reproduces the meson wave equation. We have ended up 
with an open bosonic string that moves in one dimension but whose center of mass is 


4 
















stuck. The nonlocality of (CD reflects the integration over the stringy modes labelled by 
n and the equations resemble those of There likely are other ways to localize the 

eigenvalue equation, and it is unclear whether these various localizations would be in some 
sense geometrically equivalent. 

The non-locality we have analyzed seems intrinsic, in the sense that no representation of 
the equation is known which would make it local without the addition of extra independent 
variables. There is one case when the non-locality can be neglected to a good approximation; 
When there is one very light particle, as a result of 7 << 1, the nonlocality amounts to 
determining a parameter in a local description of the light particle, as one would expect. 
We can approximate then o by assuming << 1 and then, in the light quark limit, one 
gets the usual form 

(7) 




The equation producing this result is a well known case of an explicitly factorizable 
Schrodinger equation, with a shape invariant potential 

+ («> 

This equation is an approximation, valid for the lightest meson, in the limit of light quark 
mass. The roles of parameter and eigenvalue of meson mass and quark mass got reversed. 
The equation can be thought of as an analogue of an AdS/CFT motivated partial differential 
equation determining the pion mass for given quark mass, but only at leading order in an 
energy scale given by the string tension. 

A systematic expansion in (this is not a chiral Lagrangian expansion) allows the 

calculation of higher order terms in If we knew by some other means the entire series 
in y/j expressing we could iteratively invert this to determine the nonlocality of 

the differential equation, once the lowest order (jHD is assumed as given. The exact function 
reproduces all the sigma model corrections to the leading form which here is given 
by a simple potential problem in one dimension. 

If we only know numerically, we are limited in the extend to which we can make 

inferences on the equation determining the meson masses. This is the situation we hnd 
ourselves in four dimensions, where we don’t have a general formula for the lowest order, 
starting point, equation either. 
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All is not lost though, and one can try something much less ambitious than determin¬ 
ing the full equation for the mesons, something which nevertheless carries some nontrivial 
information. Let us see what the approximate equations tells us about the function /i^(^y 7 ). 

That /i^(^y 7 ) vanishes linearly in as 7 —> 0 follows from general held theoretical 
considerations. Because we are at inhnite N, all chiral logarithms are suppressed. By a 
simple calculation we obtain: 

f‘^(7) = -^V^+4(l-^)7 + ... (9) 

Demanding the subleading correction to be smaller than 10% of the leading term gives 

fi-ji O 

However, the more natural expansion of the full ’t Hooft equation is in (^)^. The 
correction coming from the subleading term ^ is the ^ term in ©, which is only ten 
percent of the order 7 correction. Had we used this as a correction when estimating the 
accuracy of o there would have been no restriction, since the right hand side would have 
become unity, and this is about as high as the ratio can ever get. Thus, one would say that 
(0 is uniformly valid to ten percent accuracy. The factorization of (jHl) into —A^A, with A 
of hrst order in induces one to replace fA by a variable A: 

( 11 ) 

or. 


Next, one expands A in 

^ = 77-^7- (13) 

Both expansions require the inclusion of the same amount of extra (^)^”' terms at any hxed 
order, but the derivative expansion delivers reliable information about its accuracy while 
the series expansion does not. Taking for example 7 = 1, and comparing o with the exact 
result, the approximation is seen to be very good numerically, much better than the relative 
magnitude of the order 7 term would have indicated, but in agreeemnt with the order ^ 
correction. 


% 2^/3 
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TT 


71 
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The quadratic relation between the meson mass and the A variable is reminiscent of 
calculations done in the AdS/CFT context, but the meaning of A differs from the conformal 
case. In the AdS/CFT case, the quadratic nature of the relationship reflects its origin from 
an equation that has no derivatives higher than second. This local structure of the equation 
is a direct expression of the relevance of the extra (fifth for four dimensional space-time) 
coordinate of AdS. 

The approximation in which we introduce A to replace /i^ and set A = includes the 
correct leading asymptotic behavior of the quark mass dependence of the pion both at very 
small and very large quark masses. However, when we start expanding A in a power series 
in y/ 7 , we are only reproducing correctly subleading terms at the low quark mass end. Using 
the exact equation in the nonlocal form given in (HD, it is possible to set up an expansion of 
in 7 which is valid as both are very large. To this end one needs to carry out a canonical 
change of variables, from (s, to (^, —g). After that one can scale variables so that the 
heavy quark limit is smooth and one obtains a leading form of the wave equation that has a 
potential proportional to |g| and a kinetic energy of normal form. The equation describes a 
point particle in a linear potential and the fact that this holds in q, a variable conjugate to 
s, perhaps makes it easier to accept the previous picture of an open string with Neumann 
boundary conditions as describing a meson made out of massive quarks. From the new 
equation one extracts the following fact: 

,2 /;x 7 r\ 2/3 


lim 

7—^00 
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a 


= 0 


(14) 


a is the first zero of the derivative of the Airy function, a ~ —1.0188. As a result 


III, 






2 ) 
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(15) 


a behavior that cannot be represented by a truncated power expansion of 7 in A. The 
expansion in the derivative terms is inappropriate at rUg —>■ cxo. At heavy quark 

masses, a different expansion, which produces the above result at leading order, can be used 
to generate a series of subleading terms. 

From this discussion we extract the message that a parametrization of the pion mass in 
terms of a quadratically related A produces an expansion of the functional dependence of the 
pion mass on light quark mass which converges better. While the chiral expansion controls 
the structure at small quark masses, what happens in the regime of intermediate quark 
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masses already depends on the specific dynamics of the model. Because of he dynamics of 
this specific model, as opposed to other models with exactly the same chiral Lagrangian to 
order a quadartic formula for in terms of rriq has a larger domain of applicability 
than can be argued on the basis of the chiral Lagrangian alone. 

Finally, this is something we can look for by numerical means. It is not very restrictive, 
nor very solidly formulated, but testable in planar QCD in four space-time dimensions. 


NUMERICAL COMPUTATION OF THE PION MASS 


Recent work [1^ has established that planar QCD on an Euclidean torus of size and 
shape is /-independent so long as I > Ic- This is a string-like property. One refers 
to this property as “continuum reduction”, on account of the elimination of the infinite 
volume factor from the total number of degrees of freedom. For a fixed bare gauge coupling, 
g^N = Ao, this means that computations in the large N limit of four dimensional QCD can 
be performed on an L'^ lattice of relatively small size, with L > Lc{b), where b = ^. It is 
sufficient to pick L just slightly above Lc{b), since continuum reduction implies that there 
are no finite volume effects, once N is large enough. One should think about L as setting 
the minimal length scale in the problem, in this case given by 

The lattice gauge field configuration consists of a collection of SU (iV) matrices, associated 
with a link in the direction /i emanating from a site x and denoted by U^{x). These matrices 
are generated with a probability given by Wilson’s plaquette action, with coupling b, and the 
configuration is probed by a lattice fermion propagator that has exact chiral symmetry at 


zero quark mass, known as the “overlap” m propagator. Given the lattice gauge field Un{x) 
on an lattice at some coupling b, the lattice quark propagator using overlap fermions is 
denoted by G{Ufj,,mo) where 


G{U^,mo) = - - 

1 — m, 

and nio is the bare overlap quark mass parameter. is the Wilson lattice Dirac operator at 
mass mo with a so called r-parameter set to unity. The meson momentum is implemented by 
changing the gauge fields felt by the constituent quarks by independent U{1) phase factors; 
this is the so called “quenched momentum prescription”. Meson propagators are computed 


l + mo + {l- moh^^HwiUi,)] 


- 1 


(16) 







using this quenched momentum prescription and are given by 


Mvip, mo) 


Tr 


STG{Uf,e 


iPfx 

2 


,mo)Sr^G{U^e 



( 17 ) 


The definitions of S and T follow below: We choose T = 75 for a pseudoscalar meson and 
r = 1 for a scalar meson. The two quark propagators see gauge fields that differ by a 
U{1) phase and this difference is the momentum that is carried by the meson. The meson 
momentum is taken along a lattice axis: 


Pm = 


(18) 


0 if /i = 1, • • •, 3 
^ if/i = 4; 0<n<iV 

S smears the operator in the remaining, perpendicular, = 1, ■ ■ ■ ,d — 1 directions. This 
“smearing” creates an extended object that more closely approximates the true constituent 
quark structure of the meson. The smear operator is chosen to be the inverse of the gauged 
laplacian. 


1 


d-l 


s’‘ = 5E(2-7;-7), 


(19) 


M =1 


where is the gauge covariant translation operator defined by the action (T^'ip){x) = 
U^{x)'ip{x + fi). Because of fluctuations, S actually has a sizable gap and should be thought 
of as a lattice version of the covariant Laplacian with a mass of the order of the inverse lattice 
spacing. The lattice theory is in the confining, /-independent phase, when L > Lc{b)] there 
the Zjf symmetries associated with the Polyakov loops in the four directions are unbroken. 
Along with the gauge symmetry, one can use this to show 


Mr{p, rrio) = Tr 


,mo)ST^G{U^G‘^^-'^ ,mo) 


( 20 ) 



( 21 ) 

if /i = 4; 0 < n < iV, 

making it explicit that the meson propagator depends only on the difference between the 
two phases seen by the two valence quarks. 

Both the pseudoscalar meson and the scalar meson propagator computed above are given 
by sums over an infinite number of poles. Smearing reduces the residues associated with the 
higher poles. We can further reduce their contribution by working with the sum. 


1 

F{p,mo) = - [Mi{p,mo) + My^{p,mo)] 


( 22 ) 
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Excited pseudoscalar mesons and scalar mesons get closer in mass as one looks at high 


excited levels in QCD and one gets cancellations between the higher states contributing 
to F(p, mo). In practice, a single pole £t to the gauge held average of F(p, mg) works well. 

rl{mo) 


{F{p,mo)) = 


(23) 


p2 rnKrUo) 

We compute the meson propagator using a stochastic estimate of the trace in (IT7j) . We 
start by picking one Gaussian chiral source, \q >, such that 'j 5 \q >= \q >■ Then we compute 
SG{e~,mo)\q > and G{e~~,mo)S\q >. The action of the overlap propagator on a vector 
is computed using a standard multiple mass conjugate gradient algorithm. The hermitian 
Wilson-Dirac operator, has a substantial spectral gap for large N gauge fields due to 
the gap in the eigenvalue distribution of the single plaquette parallel transporter. This 
eliminates the need to project out approximate zero modes from H^, making the calculation 
straightforward. A satisfactorily accurate action of e(hf^) is achieved by just using the 21®* 
order Zolotarev approximation. The stochastic estimate of F{p,mo) is then given by 


F{p, nio) =< q\G\e , mo)S ^ \ G{e , mo)S\q > 


(24) 


A multiple mass conjugate gradient algorithm for the calculation of G easily traces the 
dependence on quark mass. Statistical errors are reduced by using the same Gaussian chiral 
source all all momenta. This numerical procedure results in an estimate for (F(p, rrio)) where 
the values at all p and rrio are correlated. 

The stochastic estimate of F{p, irto) is then fitted as a function of p at a fixed quark mass 
to give m^(mo). Minimization hts are performed using the full correlation matrix and errors 
are estimated using the jackknife method. We obtain in this way the pion mass as a function 
of the bare overlap quark mass parameter at various gauge couplings. This result is free of 
hnite volume effects since we use large enough N, leaving only hnite lattice spacing effects 
to worry about. 

In order to take the continuum limit, we hrst need to convert the bare overlap quark mass 
parameter to a physical quantity. As is well known, the quark mass itself can be defined to 
renormalize multiplicatively, and with the help of the overlap propagator this can be done 
also on the lattice. To eliminate this renormalization effect we consider instead the quantity 
moS(6) where S(6) is the bare chiral condensate at a coupling of b and zero quark mass. 
We can then convert both the pion mass, m.^-, and mo'F{h) into dimensionless quantities 


10 




b 

L 

N 

Lc{b) 


0.345 

8 

19 

4.77 

0.1675 

0.350 

8 

23 

5.97 

0.1420 

0.355 

10 

19 

6.96 

0.1265 

0.360 

11 

17 

8.01 

0.1130 


TABLE I: Simulation parameters, critical box size and bare chiral condensate 


by forming m.,^Lc{h) and moTi{h)L\{h). Plots of the functions relating these dimensionless 
quantities obtained at different lattice spacings should approximately fall on a common curve 
that ceases to depend on the lattice spacing as one approaches the continuum limit. 


NUMERICAL RESULTS 


We work with four different couplings as shown in Table d The associated values for L 
and N used in the numerical simulation are shown in the same table. The critical sizes, 
Lc{b), are known from our previous work on continuum reduction and is given by 


= e{b) = ^{TrU,Jx)) L,(6) = 0.26 


121 2i-iT^b, 


e 11 


(26) 


SO have a good estimate of the bare chiral condensate at 6 = 0.350 from previous 
1^ . but the estimates at other couplings, like b = 0.355, are not as accurate. Therefore, 


We a' 
work 

we use the known value at 6 = 0.350 and adjust the values at the other couplings such that 
all plots of m^L‘l{b) as a function of moT,{b)L‘^{b), for different values of b, fall on one curve. 
This did not have to happen, but it does, indicating that lattice spacing effects are below 
our statistical errors. The values of S(6) so obtained are shown in Table HI The value at 

n 

b = 0.355 in Table H] is higher by less than 5% compared to the numbers in but this 
is well within the error. The value at 6 = 0.345 used here is compatible with the value 
found at 5 = 0.346 in ll^ The value at b = 0.360 used here indicates that the value at 


b = 0.3585 found in 


15j | probably could go up a little bit if one does a careful analysis. 


There are potential order differences between the condensate values listed in Table H] and 
those determined using chRMT in The running of S(6) is given up to one-loop by Q 

E(6)L(!(6 )=Clln(M6)A)|* (26) 
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m„i;(b)L‘‘(b) 


FIG. 1: A plot of pion mass as a function of quark mass in dimensionless units. 

We find that C = 0.828 and A = 0.268 fits the numbers in Table HI quite well for b ^ 0.345 
and we hud a 8% deviation at & = 0.345. 

The plot of rn^Ll{b) as a function of mo^{b)L'^{b) is shown in Fig. ^ The data is then 
htted to 

A = i|v/mJL?(t) + A;-AJ; = A(A + A,) (27) 
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with 


A = m,S{b)Lt{b) + -^mll.^(b)Ll{b) + .... 


(28) 


= 6.91 and Aq = 1.03 fit the data over the whole range as shown by the solid black cnrve 
in Fig. n A plot of a trnncated chiral expansion, 



(29) 


shows that the second term in the above eqnation makes a large contribntion for 


moS(6)L^(6) > 0.2. Keeping two orders in the chiral expansion also does not agree with the 


data as well as the A parametrization. One shonld note that A >> A^ in the large qnark 
mass limit and this is not the case for the range plotted in Fig. ^ 

The parameter A^^ governs the low qnark mass regime throngh the relation ^ 


Using, 1/lc = Tc = 264 MeV, we get = 71 MeV. This translates to f,, = 123 MeV 
for SU(3) at zero qnark mass, significantly higher than the conventional valne of 86 MeV. 


Apparently, ^ effects on the pion decay constant are larger than on the glneball mass or the 


hnite temperatnre phase transition, all measnred in nnits of string tension. 

SUMMARY 

For simple qnantnm mechanical problems the spectrnm and scattering data can be nsed 
to determine the potential by the inverse scattering method. Conld we nse spectral meson 
data at inhnite N to determine the string descriptions of mesons ? 

In this paper we took a very primitive first step with this philosophy in mind. We were 
led to the specnlation that the strnctnre of QCD in the planar limit, where chiral logarithms 
are snppressed, explains why mass formnlae of the type = 2Bmq + cm^ seem to work 
even when is not small. Of conrse, this may be just saying that the heavy quark regime, 
where mjr ~ 2mq is smoothly connected to the light quark regime. We observed that this 
could also reflect the approximate validity of an equation similar to (0. 
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